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Abstract 

We discuss different aspects of the Higgs self-interaction in the MS and the on-mass- 
shell (OMS) scheme. The running coupling A(/x) is investigated in great detail. The 
three-loop coefficient of the /3-function in the OMS scheme is derived, and the three- 
loop running coupling is calculated. The breakdown of perturbation theory for large 
Higgs masses Mh is analyzed in three physical observables for which two-loop re- 
sults are known. Requiring the dependence on the renormalization scale to diminish 
order-by-order in A, we find that perturbation theory breaks down for Mh = 0(700 
GeV) in Higgs decays. Similarly, Mh must be smaller than 0(400 GeV) for pertur- 
batively calculated cross sections to be trustworthy up to cm energies of 0(2 TeV). If 
the Higgs sector shall be perturbative up to the GUT scale, the Higgs must be lighter 
than 0(150 GeV). For the two-loop observables examined, the apparent convergence 
of the perturbation series is better in the OMS scheme than in the MS scheme. 
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I. INTRODUCTION 



One of the least tested sectors of the Standard Model of elementary particle physics is the 
Higgs sector generating the masses of all particles via the mechanism of spontaneous symmetry 
breaking. This mechanism implies that the self-coupling A of the Higgs particle is proportional to 
the square of its mass, Mh- Hence a heavy Higgs particle may cause the breakdown of perturbation 
theory in A. 

When calculating decay rates or cross sections beyond the tree level term of the perturbation 
series one must specify the renormalization scheme to define the coupling constants and the par- 
ticle masses appearing in the analytic expressions. In the Higgs sector one usually adopts the 
on-mass-shell scheme (OMS): Here the squared mass coincides with the physical pole of the prop- 
agator, the vacuum expectation value v = 246 GeV of the Higgs field is renormalized such as to 
cancel tadpole contributions, and the coupling is chosen to satisfy 

AoMs - ^ - (1-1) 

to all orders in perturbation theory Here Gi;' is the Fermi constant. When discussing the per- 
turbation series of physical observables one must first distinguish processes involving two largely 
separated mass scales from those containing only one scale, Mh. Prototypes of the first species are 
cross sections at LHC energies. They contain potentially large logarithmic terms Xln^y/s/Mn) 
with y/s being the energy of the process. These terms may spoil the smallness of radiative correc- 
tions. Yet they can be summed to all orders in perturbation theory with the help of renormalization 
group (RG) methods. The corresponding equation in the OMS is the Callan-Symanzik equation, 
which can be solved in the limit of large -y/i when all particle masses can be neglected. The second 
kind of physical observables are one-scale processes. Examples are two-body decay rates of the 
Higgs particle into (almost) massless particles. We will see in the following that RG methods are 
also a useful tool to judge the accuracy of perturbative results in these cases. 

Due to dl . 1| ) the use of the OMS appears natural, because Aqms is directly related to physical 
observables. Yet it is also useful to consider mass independent renormalization schemes, the most 
prominent example being the MS-scheme. The reasons are the following: 

i) Mass independent schemes allow for an exact solution of the RG equations, i.e. mass effects 
can be systematically included. 

ii) The analysis of scheme dependences provides a test of the reliability of perturbation theory 
since the results obtained to order A" in different schemes formally differ by terms of order 

A"+i. 

iii) Results obtained in mass independent schemes involve an arbitrary parameter, the renormal- 
ization scale jJL. For perturbation theory to work it is necessary that the dependence of the 
result on fi diminishes order by order in perturbation theory. We will use this fact extensively 
in the discussion of the breakdown of perturbation theory. 



^We consistently neglect all gauge coupling corrections. Ar is hence taken to be zero. 
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The paper is organized as follows: In the following section we discuss the running coupling con- 



stant in the two-loop approximation. In Sect. we give the two-loop relation between the quartic 
Higgs coupling in the OMS and the MS scheme. We calculate the three-loop OMS coefficient of 
the [3 function as well as the coefficients of the leading and next-to-leading logarithms. Subse- 
quently, we investigate the scheme and scale dependence of two-loop physical observables using 
the three-loop running coupling. Special attention is paid to the breakdown of perturbation theory 
for heavy Higgs masses. We start by examining the bosonic and fermionic Higgs decay widths 
(Sect. ^|). Finally, we look at scattering processes which involve the longitudinally polarized 
gauge bosons and Z^. These scattering processes give the most stringent bounds on a pertur- 
bative Higgs mass (Sect. 0). 



II. PERTURBATIVE RUNNING COUPLING AT DIFFERENT ORDERS 

We first discuss the coupling constant A(/u), whose running is encoded in the /3-function. To 
three loops, the beta function of the Higgs quartic coupling is defined as 

^ ^d^i 167r2 (167r2)2 (leTr^)^ \ ) ^ ' 

We neglect all contributions from gauge and Yukawa couplings. This is an excellent approximation 
for large values of A. To two loops the coefficients are [[T|] 

/?o = 24, A = -13/?o = -312. (2.2) 

The three-loop coefficient /?2 is scheme dependent. We restrict the discussion of this section to the 
scheme independent two-loop results and return to (32 in Sect. |n[ 

Equation ( |2.1| ) is valid in any mass independent scheme with fi being the renormalization 
scale accompanying the coupling constant in the Lagrangian. In the OMS scheme, the Callan- 
Symanzik equation describes the response of some Green's function to the scaling of its external 
momenta according iopi /i//xo ■ Pi which is related to the corresponding scaling of the coupling 
given by Eq. (pTTj). The values of /5o and (3i are scheme independent, so that we don't choose 
a specific renormalization scheme until later. The determination of A(/i) proceeds in two steps: 
First at some initial scale fiQ ~ Mr the coupling A(/io) is obtained from ( |LT| ) or an equivalent 
relation, if the scheme under consideration is not the OMS scheme.0 Then ( j2.lD is solved for 
= \[\(fio), fi/ Ho]. The question whether perturbation theory works in a specific physical 
process therefore depends on the two parameters A(;Uo) and h/ Hq, which are related to the phys- 
ical parameters Mh and the energy of the process. The cause of a possible breakdown of 
perturbation theory is twofold: First, the larger the Higgs mass, M^, the larger is A(;Uo) due to 



''For the MS scheme, the expression for A(/io) is derived in section see Eq. ( |3.3[ ). 
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dlTTI). Second, the resummation of possibly large logarithms, \n(^/s/MH), introduces the running 
coupling A(/i) with fi ^ y^, and X{fj,) increases for increasing /i ^ y/s. 

Defining A at some initial scale ^q, the solution of ( |2.1D yields A(/i) at any other scale yU. The 
expansion of X(fi) around A(/io) shows that A(/i) resums powers of the logarithm ln(/i//io) mul- 
tiplied by powers of A(yUo)- The coefficients of the leading logarithms (LL), A"^^(/io) ln"(/i//io)> 
depend only on (3q, those of the next- to-leading logarithms (NLL) depend on both Pi and Pq, and 
so on. We will give these coefficients explicitly later. The first point we want to stress is that 
beyond the one-loop approximation there are several different solutions of Eq. ( p] ) due to the 
truncation of the perturbation series. All of them agree to the order at which the perturbative series 
of the beta function, (]2.1D, is truncated, but they differ by terms of the neglected order. For a viable 
perturbative treatment, these higher-order differences should be negligible. We now look at four 
different solutions for the running coupling A(/i). 

Setting simply the coefficients of the neglected terms in ( |2.1D equal to zero and integrating the 
remaining expressing exactly, we obtain an implicit equation for A(/i): 

^ = 1 - /5oA(/.o) In f I A(/.o) In ( (AlIAA^M/^) (2.3) 



A(/i) ' Uo/ Po '"^'^ \{Po + PiX{f,o))X{fi)) 

in the NLL approximation, and 



PoPlXifio) 



Xb — Xa 



J_ In ~ '^^^^ - — In ~ 

Aft Xb - A(yUo) A2 Aa - A(/io) 



(2.4) 



with 



Xa/b=-^±yp!-w2 



in the next-to-next-to-leading log (NNLL) approximation]^ Here ( |23| ) is the two-loop result pub- 
lished previously [^. In ( P3| ) and ( jZ^ ) the useful abbrevation A = A/(167r^) has been used. To 



obtain A(/i) one has to solve ( p3| ) or ( |2.4D by numerical methods. Direct numerical integration of 
the original equation, ([2.ip, yields the same result for A(/i). Let us call this form the naive solution. 
It would be the exact result if the neglected coefficients of the /3-function were really identical to 
zero. Yet ( |2.3D and ( [2.4[ ) contain logarithmic terms belonging to the neglected higher orders of 
RG improved perturbation theory. E.g., the naive NLL result for X{fi), ([2. 3D, partially resums 
NNLL logarithms A""'"^(/io) ln"(/i//io)- However, there are further NNLL terms from irreducible 
three-loop contributions. Therefore, the naive NLL result contains an inconsistent resummation 
of NNLL terms. Similarly, Eq. ( |Z4l ) includes an inconsistent resummation of higher-order terms, 
even though it is correct to the order considered. This fact is true whenever beta functions are 
integrated numerically at and beyond two loops. 

For a consistent result, one needs to ignore terms of the neglected order in A when integrating 



(2.1). This yields the consistent NNLL result 



^Setting (32 = 0, Eq. (2.4) reduces to (2.3), of course. 
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A(/io) 
A(/x) 



(2.5) 



The consistent NLL result is obtained from ( [2.5[ ) by dropping the underlined terms. The final result 
for the consistent solution is obtained by solving Eq. ( |2.5| ) numerically. 



Third we solve ( [2.5D iteratively by first substituting the one-loop result for A(/i) into the RHS 
of the equation above and then repeating this step with the result of the first substitution. This 
yields the iterative answer 



A(/x) = A(/xo) < 



l-/3oA(/io)ln 



+ ^M/^o) In 
Po 



1 - /5oA(/io) In — + ^A(/io) In ( 1 - /3oA(/io) In — 
/^o Po V fJ'O, 



Po /io V Ato, 



-1 



(2.6) 



We stress that no further expansions in A(/io) are possible, because each A(/io) multiplies a large 
logarithm ln(/i//io). Eq. ( |2.6D is the NNLL result, and the NLL expression is again obtained by 
dropping the underlined terms. 

The fourth solution of ( [2. ID is constructed in analogy to QCD: The integration constant ob- 
tained in the integration of ( [2.1[ ) can be absorbed into a scale parameter Ah'- 



167r2 



PoHf] 



2^1n[ln(^)] 4/3f Inlln 



+ 



PI 



+ 4 



P0P2 - PI 



Pt 



-4 



O 



ln=^[ln 



(2.7) 



We call this result the QCD-like solution. The above equation defines the scale parameter Ah, 
and it is written such that no term of the form const./ ln(A^//i^) appears in the square brackets. 
This is identical to the definition of the QCD scale parameter Aqq^ [||]. The definition of the NLL 
parameter A^^^ is obtained by dropping the underlined terms. The above definition of Ah works 
in any renormalization scheme. The actual numerical evaluation of Aj^ depends on the boundary 
value A(/io)5 which is scheme dependent. If for example A(/i) on the LHS of ( l277| ) is given in the 
MS-scheme, Ah on the RHS equals Af^. 

There are two important differences between the Higgs sector and QCD: First, Eq. ( |277| ) holds 



for fi <C Ah, while in an asymptotically free theory like QCD the analogue of ( p77| ) is valid 
for yu ^ Aqcd. Second, Aqcd is the only fundamental scale parameter of QCD (with massless 
quarks), while in the Higgs sector A// is related to the Higgs mass (see Table |). Yet if the Higgs 
mass turns out to be large, one will have to parametrize non-perturbative effects in terms of an 
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effective Lagrangian, and Ah will be the natural scale entering the effective couplings. In QCD 
perturbation theory breaks down for < 3 Aq^^ and binding energies equal a few times Aqcd- In 
the Higgs sector the breakdown of perturbation theory likewise occurs for fi > A///3. In Table | 
we tabulate the relation between Mh and A^. When the OMS scheme is adopted, in Table | 
is the physical mass (see (|l.ip). In other schemes the tabulated values for Mh correspond to the 
tree-level relation \{fio = M^ff^ ) = {M^" )^/(2f ^) and — in contrast to the OMS scheme — 
radiative corrections to this relation have to be taken into account in order to obtain the physical 
Higgs mass from M^^^ . For the MS-scheme this relation is encoded in ( |3.3| ) of Sect. |II[ 





A^^(GeV) 


A^^^^(GeV) 




AS(GeV) 


MZ."(GeV) 


100. 


4 ■ 10^6 


6 ■ 1037 


7 ■ 1037 


7 ■ 1037 


2.2 - 1036 


200. 


9 ■ lO^o 


6 • 10^1 


7-10" 


7-1011 


5.1 - 10^° 


300. 


2.1 ■ 10^ 


8.7-106 


1.0- 10^ 


1.0 • 107 


1.2 - 10^ 


400. 


5.8-10^ 


1.8 ■ 10^ 


2.0 ■ 10^ 


2.1 • 10^ 


3.4- 10^ 


500. 


1.2 ■ 10^ 


2.9 • 10^ 


3.3-10^ 


3.6 ■ 10^ 


7.0 - 103 


600. 


5.51 ■ 10^ 


1.13 ■ 10^ 


1.25 ■ 10^ 


1.40 ■ 10^ 


3.2 - 103 


700. 


3.57-103 


6.45 ■ 103 


6.92 ■ 103 


7.99 ■ 103 


2.1 - 103 


800. 


2.78 ■ 103 


4.52 ■ 103 


4.71 ■ 103 


5.64 ■ 103 


1.6 - 103 


900. 


2.41 ■ 103 


3.57-103 


3.60- 103 


4.52 ■ 103 


1.40 - 103 


1000. 


2.22 ■ 103 


3.03 ■ 103 


2.94- 103 


3.91 • 103 


1.29 - 103 


1100. 


2.13 ■ 103 


2.68 ■ 103 


2.48 ■ 103 


3.57-103 


1.24- 103 


1200. 


2.09 ■ 103 


2.41 ■ 103 


2.18 ■ 103 


3.37-103 


1.21 - 103 



TABLE I. The values for the scale parameter Ah obtained from ( ^?7| ) for 
left column corresponds to M^ff^ = 



M]f in the 



2A(/io 



For A 



AoMS this relation receives no radiative 
corrections. If the coupling is defined in the MS-scheme, use (p^) first to calculate Aj^(^o 



ALL 



from the physical Higgs mass to the desired order. Then use the previous equation to obtain 
corresponds to the location of the one-loop Landau pole. The difference between the fourth and fifth column 
is caused by the scheme dependence of j32- The last column contains the two-loop values of ix for which the 



iterative solution, (2.6), assumes its maximum. 



We now completed the definition of the four different solutions for the running Higgs coupling: 
naive, consistent, iterative, and QCD-like solution. It is interesting to note that at one loop all four 
solutions are identical.^ This is clearly not the case anymore at two-loop and higher orders, and 
we will discuss the differences below. Let us stress again that all solutions correctly sum the large 
logarithm ln(/i//io) within the calculated order. The difference between these solutions is of the 
neglected order as can be seen when expanding the solutions in A(/io)- If perturbation theory is 
applicable this difference should be numerically small, giving a very simple criterion to find the 
values of Mh and ji beyond which perturbation theory clearly fails. 

In Fig. H], we compare the fi dependence of the four solutions at one and two loops. We 
choose the three mass values Mh = 200, 500, and 800 GeV. To obtain a meaningful compari- 



^To see this, set /3i = /32 = in the four different solutions. 
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son between the one- and two-loop results for a given value of Mh, we take for both orders the 
same expansion parameter Aq = A(;Uo), and choose it to be defined by the tree-level relationship 
Ao = {M^'^'^)'^ /(2f^). This has the additional advantage that we do not yet have to specify a 
renormalization scheme in Fig. |T]. Since the first two coefficients of the /5-function are scheme 
independent, the whole scheme dependence resides in the relation between the tree-level mass 
(which labels the different curves in Fig. P and the physical Higgs mass. In the OMS, the two 
mass definitions coincide. 
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FIG. 1. Different solutions of the one-loop and two-loop RGB equations for the Higgs quartic coupling 



A: the one-loop solution (long dashes), the two-loop naive solution according to (2.3) (long-short dashes). 



the two-loop consistent solution of ( |2.5[ ) (solid line), the two-loop iterative solution ( p.q ) (short dashes), and 
the two-loop QCD-like solution expanded in powers of 1/ ln(A|^//i^) given in (2.7) (dots). The crosses 
show the normalization point: A(^ = Mf^^^) = Aq- To allow for a meaningful comparison of the different 
orders in perturbation theory, we choose Aq to be the same at one and two loops. The curves are labeled 
by the tree-level Higgs mass (in GeV), corresponding to the tree level relation Aq = {JVPff'^)'^ / {2v'^). In 
renormalization schemes different from the OMS, radiative corrections to this relation must be included 
(see text). 



At one loop, the four solutions are identical (long dashed curve). The coupling X^^{fi) ap- 
proaches infinity as /i ^ A^^, which is usually referred to as the Landau pole [0. At two loops 
the breakdown of perturbation theory is clearly visible in the two-loop curves for Mh = 800 GeV. 
It manifests itself in a very different behaviour of the four solutions, although they all are defined 
to have the same value for jiQ = Mh- The one-loop Landau pole has completely vanished from all 
two-loop solutions. The naive (long-short dashes in Fig. |T]) and consistent (solid) solutions always 
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have positive slopes in accordance with the positiveness of the perturbative /^-function, Eq. (|2.1|). 
The consistent solution can be looked at as the "average" of the one-loop solution and the two- 
loop naive solution. The naive solution approaches the fix-point value Aj^^^ = 12.147, the zero of 
the two-loop beta function. This simply reflects the fact that one cannot use perturbation theory 
to gain information about large values of the coupling. This is even more evident from the other 
two solutions. The QCD-like result for has a global maximum. Interestingly enough, 

for values of Mh > 1263 GeV, the value Aq = Mjj/{2v'^) is larger than the maximal possible 
value of the QCD-like running coupling, so that X{n = Mh) = Aq can no longer be satisfied. The 
unphysical maximum of the QCD-like solution appears for /i^^ = and perturbation 

breaks down well before is equal to Kh- The iterative solution features a maximum at 

which is located at 1615 (7033) GeV for Mh = 800 (500) GeV. The two-loop iterative solution 
approaches zero at the same value of /i at which the one-loop solution has its Landau pole. 

The fact that some perturbative solutions of the two-loop running coupling show a maximum 
gives a first criterion for the breakdown of the two-loop perturbative treatment. From the definition 
of the P function, Eq. (|2.ip, we know that in the perturbative regime the slope of the running 
coupling has to be positive as the one-loop term dominates over the two-loop term. If the iterative 
solution for A(/i) has negative slope, then perturbation theory is not valid anymore. Hence Eq. ( |2.8| ) 
is possibly a measure of the range of /i for which perturbation theory is meaningful. 

The previous discussion of the four solutions is valid for any value of Mh- Yet if we restrict 
ourselves to values /i < 5000 GeV, features like poles and maxima are only apparent for Mh = 
800 GeV, and not for 500 and 200 GeV (see Fig. [1]). For the lower mass values one has to go to 
(much) larger values of /i to observe the breakdown of perturbation theory. In fact, if Mh = 200 
GeV all four solutions of the running couping agree extremely well for n < 5000 GeV. Choosing 
Mh = 500, however, already results in an uncertainty between the different two-loop solutions 
of about 25% at = 5 TeV, and the one-loop solution is more than 50% larger than the two-loop 
solution at /i = 5 TeV. 

If the Standard Model Higgs sector is to remain perturbative up to much higher energy scales 
the restrictions on Mh get much more severe [^. To avoid the Landau pole and keeping the 
coupling A perturbative, one can define an embedding scale below the Landau pole and require 
the running coupling to be smaller than a certain value. For example, taking the embedding scale 
to be 10^° GeV the authors of implement the "automatic fixing procedure" requiring (in our 
notation) A(10^°GeV) < tt^ and find an upper bound of Mh ~ 230 GeV when neglecting the 
top-quark Yukawa coupling. Using Eq. (Ob and requiring /zj^^gjj^ to be equal or larger than 10^° 



GeV, we find the upper limit Mh ^ 210 GeV. 

The analysis of the running coupling alone can clearly only yield a necessary criterion for the 
validity of perturbation theory. Taking Mh = 800 GeV, we can be sure that perturbation theory 
breaks down for any cross section with i/i ^ 2 TeV, because this energy is too close to the scale 
parameter Ah and already larger than the corresponding maximum of yU, ix^^^-^{Mh = 800GeV) ~ 
1600 GeV. For Mh < 500 GeV, the perturbative RG treatment of the running coupling could 
be trustworthy up to values of /imix ~ ^ TeV, sufficient for studying Higgs physics at the LHC. 
Nevertheless, choosing Mh = 500 GeV the running coupling X^^^{fi = 7TeV) = 4.4 is sizable. 
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To make a final judgement on whether perturbation theory works for Mh = 500 GeV and 
TeV-energies, one must in addition investigate the perturbation series of the physical process of 
interest. As we will see later on, the perturbative solution of the running coupling may show a 
reasonable convergence, but the numerical value of the running coupling is already too large to 
calculate physical quantities like cross sections and decay widths in perturbation theory. This will 
be the subject of sections |^ and M. 



III. SCHEME DEPENDENCE: OMS VS. MS FORMULATION 



So far the discussion of the running coupling has been independent of a special renormalization 
scheme. We have already listed the reasons for studying different renormalization schemes in the 
introduction. 

In the following we will look at the MS and OMS scheme, and examine the scheme dependence 



of the coupling and the three-loop coefficient P2 in (|2.1D. Before going into detail we would 
like to remark two points: First, the OMS coupling Aqms is to all orders in A directly related to 
measurable quantities (the muon lifetime and the Higgs mass) via (|1.1|). This is not so for any other 
renormalization scheme, where the RHS of ( |1.1[ ) receives radiative corrections, which depend on 
an additional parameter, the renormalization point Hq. Second the scheme dependence of P2 makes 
the coupling run differently in different schemes. One can in general adjust the scheme such as to 
achieve any desired values for /5„ with n > 2. A criterion for a "good scheme", however, cannot 
be founded on the smallness of the running coupling alone. Instead one has to consider physical 
observables, in which the coefficients of the perturbation series depend on the scheme as well. 
This will be done in the following sections. 

The starting point of the analysis is the two-loop relation between the bare and the renormalized 
coupling: 



Pol 
2 e 



+ ci + ecu + e^ci2 + O(e^) 



f ^+(/5oCi + |)^ + C2 + ec2i + 0(e2) 



+ (-Y^/^oCll + g/^oC? + -P0C2 + —PlCl + -(32) - + C3 + 0( 



. (3.1) 



In the MS -scheme the various quantities are given as 

Xmif^) 



X 



167r2 



^ = Aire 



Ci = C2 = C3 = Cii = C12 = C21 



0. 



In the OMS scheme, they instead read 



A 



VOMS 



Ml 



Cl 



25 - 371^3 = 8.676, C2 = 378.5 , Cii = 3.821, (3.2) 



and C3, C12, and C21 are unknown. Knowing the bare coupling to two loops in both the OMS 
scheme and the MS scheme, we can express the MS coupling in terms of the OMS coupling, and 
we can calculate the difference between and (3^'^. 
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A. The MS-coupling and the renormalization point /uq 



Since the OMS coupling is a function of Mh (or v) and is fixed to all orders by ([TTT|), the 
definition of in terms of Aqms is equivalent to a relation between A^jg and Mh- With Cj and 



Cjj refering to the OMS quantities in ([3 .21) we find: 



OMS 



+ 



+ 





2 

) In 








3 

1 In' 


2 > 





2 A'-Q 



+ 



3 A'-o _| 

2 



V 2 

/3o /5 



+ /?oCi In —5- + C2 - /^ocii A 



Mi 



OMS 



Pi + ^/3oCi) In 
+ ((3C2 - SPoCn) y + Aci + ^[3^'^^) In 



2 A'-Q 

Mi 



Mi 



+ C3 + ( t/^oCi2 - cicii - C21 ) l3o - t/5iCii Aqms 



3 



+ O (A^ms) + 0(e) 



(3.3) 



A 



OMS 



12 In 



+ 144 In 



Mi 

2 A^o 
Mi 



+ 8.676 AoMS 



52.22 In ^ + 286.8 j A^ms 



+ ( 1728 In^ Sr - 932.0 In^ + 9736.8 In ^ + ^o' 



Mi 



Mi 



Mi 



'^OMS 



+0(e) 



(3.4) 



which agrees to one loop with the result of Sirlin and Zucchini [|T0|]. The three-loop constant term 
d^Q depends on the yet unknown OMS coefficients C12, C21, and C3. 
For convenience we give the inverse formula of ( |J!4l ) as well: 

,2 \ 



A 



OMS 



•^Ms(/^) 



-12 In 



Mi 



8.676 AMs(/i) 



+ f 144 ln=^^ + 364.2 In 



Mi 



Mi 



136.3 Ay/i) 



Msy 



Ofel 



(3.5) 



Eq. ( p.3| ) defines the expansion parameter \{fio) of the running coupling when using the MS 
scheme. It shows that Aj^ is completely determined by specifying Mh (and thereby Aqms) and 
a renormalization point /io at which ( p.3p is imposed. The scale /io is unspecified. However, to 
ensure that the logarithms ln(/io/MH) stay small at all orders, it should be chosen of the order of 
Mh. 
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FIG. 2. The MS coupling, Aj^, as a function of the renormalization scale fi, fixing the physical Higgs 
mass at 200, 400, and 600 GeV. The value of fi at which the two-loop (one-loop) MS coupling is equal to 
the /i-independent OMS-coupling, Aoms> is indicated by a diamond (crossed diamond). 

Let us emphasize that throughout the paper Mfj denotes the pole mass, even when discussing 
MS renormalization, and the vacuum expectation value v is also chosen as defined in the OMS as in 
[|T0|]. Expressions involving the running mass of the MS-scheme can systematically be expressed 
in terms of the pole mass M^- Despite the use of the pole mass M^, the MS renormalization 
scheme maintains all the advantages of mass independent schemes mentioned in the introduction. 
For example Aj^(/^o = M^) = (800 GeV)^/(2t'^) corresponds to a physical Higgs mass of 
Mh = 720 GeV at one loop, and Mh = 681 GeV at two loops. 

In Fig. 1^ we show the MS coupling as a function of /io for different values of Mh, varying 
/io in a typical range, Mh/2 < /iq < 2Mh- We find that the MS coupling is larger than the 
OMS coupling for most of the range of fio examined. If Mh is larger than 400 GeV, the value 
of jjQ at which Aj^(/io) equals Aqms changes significantly when going from one loop to two 
loops. At one loop, the two couplings are equal if /xq ~ 0.7 Mh for any value of Mh, because 
Po/2 ■ ln(0.7^) = — Ci. At two loops, the relation gets a more complicated mass dependence. 
We find that for Mh = 400 (600) GeV the two-loop couplings of the two schemes are equal if 
/io = 0.6 Mi^ (0.4 Mh). 

At first sight it appears reasonable to choose /io in such a way that calculated radiative correc- 
tions become small or even zero. This criterion of (fastest) apparent convergence (FAC) has had 
some supporters in QCD RG analyses over a decade ago (for a criticism see [pT]]). Also in 
the one-loop FAC scale /io = 0.1 Mh, at which the one-loop correction in ( |3.3D vanishes, has been 
used to define Af^. Yet in the following we will explain why this choice of the scale is not very 
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useful, in particular in conjunction with the analysis of the breakdown of perturbation theory. We 
can already anticipate this fact from the above observation that for large values of Mh the two- 
loop FAC scale differs sizably from its one-loop value. For a more detailed investigation consider 
the generalization of the square bracket in ( |33| ) to the order Aqms" 



^Msif^o) — AoMS 



N 

n=l k=0 



dnk In 



j4_ 

Ml 



(3.6) 



i.e. one has d 
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= Ci, (ill = /5o/2 and so on. Application of the renormalization group equation 
/io) yields recursion relations for the dnk^- The solution for the coefficients dnn 



( gyiD (with /i 
of the leading logarithms is well-known 



7 _ / 




(3.7) 



Similarly we derive the general relationship 



d„ 




(3.8) 



for the next-to-leading logarithms with 



n— 1 k -1 

T.T.- 

k=ij=i J 



where ^(x) = Inr(x). For n < 3 the coefficients dnn and dn,n-i have already been given in 
(p3|), where also the NNLL coefficient d-^ involving P2 is displayed. 

We remark here that the coefficients ( |3.7| ) and ( |3.8| ) have the same form for any expansion of 
a running coupling A(/i) in terms of some ^u-independent coupling with only the coefficient ci 
changed correspondingly. For example, the expansion of A|^(/i) in terms of Aj^(M//) yields a 
series of the form ( ]J!6| ) with Aqms ~^ ^ms(^^h) and dio = ci — in (|J!8[). 

Now the choice /iq = nullifies all logarithms to all orders in (p.6|). The FAC scale hq = 



0.7 Mh instead nullifies the Aqms"'^^^'^ (|3.6|) and (Ob, but the price to be paid is the appearance 



of logarithms in higher orders. Since we know the NLL coefficients C2 and cn, we can check the 
effect of the FAC scale setting on the AoMs't^^^^ (p3|): For fiQ = Mh the coefficient of Aqms 
equals 287, but for /io = 0.7 Mh one finds the larger coefficient 324 which increases the impact 
of higher order terms, especially if the coupling is large. This explains the above finding that the 
FAC scale changes significantly for large values of Mh when passing from the one-loop to the 
two-loop order. Hence the FAC scale setting pushes large terms from the calculated orders into 
the uncalculated higher orders of the perturbation series. Yet a clever choice of /xq should yield the 
opposite and keep these higher order terms small. 

Of course one is not forced to use /io = Mh exactly. But what is the allowed range for /io? 



Clearly Mfj) in ( |3^ ) must be kept small, so that the logarithms do not become dominant in 

higher orders. A reasonable interval for /io should therefore obey 



In 



m+l 



< 



d W 



M'h 



12 



for m = 0, . . . n— 1 . Especially the size of the logarithms should not exceed the constant term | dno \ . 
With the two-loop information ( |3.4D at hand we first compare the size of both ^22 lii^(/io/^H) ^^'^ 
d2i ln(/ig/M|f) with the non-logarithmic term, ^20 = 286.8, and find that for 0.5 < fio/Mn < 2.0 



their magnitudes do not exceed \d2o\. Second, using ( p77| ) and ( p^ we can compare the leading 
logarithm dnn with the next-to-leading logarithms dn,n~i ln"~^(/Xo/M^) to all orders 

in perturbation theory obtaining the smaller range 0.8 < ^q/Mh < 1-25. At three loops we 
can use the result ( |3.1 ip to repeat the game with the next-to-next-to-leading logarithm, which is 
multiplied by rfgi = 9736.8 in (gj]): 1^33 ln^(/ig/M^)| and ln^(/io/^^^)l are smaller than 
|ci3i ln(/i2/M|r)| for 0.3 < iiq/Mh < 3.3. For ^iq/Mh = 0.3 or fio/Mn = 3.3 this term, however, 
exceeds 2 ■ 10^, which is above the value one expects for the yet unknown constant \d^o\. Since 
dss, dsi > and ^32 < the choice hq/Mh = 0.7 < 1, which has been so seductive in the order 
Aoj^g, leads to the fact that the three logarithmic terms in the order Aqms ^^d with the same sign 
to yield —8045. Moreover from ( |3.7D and ( |3.8D one realizes that for /ig < Mh the leading and 
next-to-leading logarithmic terms have the same sign to all orders Aq^jq except for n = 1 and 
n = 2! Hence the one-loop FAC choice /xq = 0-7 Mh pushes large terms into the higher orders. 
The two-loop FAC scale is even lower, increasing all higher order terms even more. 

We conclude that our lack of knowledge of the higher order terms of the perturbation series 
forces us to consider any choice for fiQ in the range 0.8 < hq/Mh < 1.25 (which may possibly 
be relaxed to 0.5 < Hq/Mh < 2.0) with equal right. Changing yUo in the n-th order perturbative 
expression for some observable changes the result by terms of the neglected order A"+^. When 
perturbation theory works the dependence on /iq diminishes order-by-order in A. We will use a 
similar criterion to find the breakdown of perturbation theory in sections |^ and |^. 

Let us close this section with a final remark on the arbitrariness of /iq: Cross sections with 
cm-energy ^/s and expressed in terms of A^(yu) involve the logarithm ln(/i^/s). Using ( |3.3| ) 
with jjL = jjLQ, this logarithm would be large. Using the running coupling Ay^(/i) evaluated at 
a scale /i ~ i/i, the logarithm is summed to all orders in perturbation theory. The arbitariness 
in the choice of /Xq ~ Mh and /i ~ reflects the fact that one can sum an arbitrary small 
constant together with the large logarithm. This feature is also present in the OMS scheme, but 
less apparent. For example, the authors of [0] sum the constant ci in (p.3|) together with ln(/i^/s). 



B. The NNLL /3-function in the OMS and MS-scheme 

Neglecting all couplings except for the Higgs coupling A, the Higgs sector of the Standard 
Model is equivalent to a spontaneously broken cj)"^ theory with N = A real scalar fields. Theories 
with spontaneous symmetry breaking have the remarkable property that the counterterms needed 
to make the theory finite are the same for the broken and unbroken symmetry, if a mass inde- 
pendent renormalization scheme is adopted [[T^]. Since the beta coefficients are calculated from 
the counterterms of the coupling, Eq. (|3.1D, the beta function is also identical for both the broken 
and the unbroken theory. Hence one can calculate the {n + l)-loop coefficient of the Higgs 
sector in two steps: First, obtain in an unbroken cj)'^ theory with = 4 real scalar fields by 
calculating the divergent parts of the four-point function to the [n + l)-loop order. This is easier 
than working in the broken theory, because only a four-point coupling is involved. Once this is 
accomplished, of the broken theory is known, too. Second, calculate the scheme dependence 
(3^^ — (3^^^- This only requires the calculation of the ra-loop (not n + 1-loop) finite parts of 
self-energy diagrams. 
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The three-loop coefficient /?2 of the (j)^ theory has been calculated in the MS scheme [ [131 ] : 

^Ms ^ 7^7Q ^ 4032C(3) = 12022.69 .... (3.9) 

To obtain Z^^'^^, we return to Eq. ( [3.1D to calculate the scheme dependence of the coefficient /?2- 
As stated above, unlike P2 itself the scheme dependent difference between and can 
entirely be obtained from two-loop quantities: 

= P^^^ - 5400 + 8688C(2) - 2160C(3) - 27367ry3 

+11527rCl(7r/3) + 5184y3Cl(7r/3) + 3888/^5 
= (3^^^ + 7784.45... (3.10) 

where = 0.92363 ... is the value of the all massive Master diagram with = Mfj, which 
has been evaluated numerically in [|]] . ( [3.100 can be obtained from the definition ( p] ) of the f3- 



function or by comparing the bare coupling expressed in terms of with its expression in terms 

of AoMS- 

Using the result for we obtain 

^OMS ^ 4238.23 ... . (3.11) 

Our analytical^ result agrees with the numerical result obtained in [ p^ ] to better than six digits. 

In Fig. we show the yU-dependence of the NNLL (three-loop) running coupling in the OMS 
scheme. The consistent solution is almost identical to its NLL (two-loop) result even for = 
800 GeV, whereas the other NNLL solutions show a behaviour very different from their NLL 
results for large Mh. For = 500 GeV and fi < 5 TeV all four NNLL solutions show a very 
nice convergence. For such a value of Mh, however, the LL running coupling is not an adequate 
approximation for the upper values of fx considered. 

In the MS scheme, the convergence of the running coupling is rather poor for Mh = 500 
GeV and above. This is also true for the consistent solution. For Mh = 750 GeV, the consistent 
NNLL solution of the MS running coupling is not defined anymore if yU > 950 GeV. The poor 
performance of the NNLL MS coupling is due to the term involving ^2- Although is only 
three times larger than the coefficient Pf — P0P2 entering the running coupling ( |2.5[ - P77[ ) is 
44 times larger in the MS scheme. This is caused by a numerical cancellation in the OMS scheme, 
where Z?^^^/ Pi ~ Pi/ Pq- We remark here that in a scheme with exact geometrical growth of the 
coefficients, Pn+i/ Pn = Pi/ Pq, the consistent solution ( p3| ) equals the consistent NLL result to 
all orders. 



IV. SCHEME AND SCALE DEPENDENCE OF HIGGS DECAYS 

The accuracy of perturbation theory and its breakdown as Mh increases can only be investi- 
gated in physical observables. Here processes in which all mass parameters are of the order of 



^apart from the numerical constant which is defined by a single Feynman diagram 
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FIG. 3. Different solutions of the three-loop RGB equations for the Higgs quartic coupling A: the 
NNLL (three-loop) naive solution according to ( |2.4[ ) (long-short dashes), the NNLL consistent solution of 



(2.5) (solid line), the NNLL iterative solution (short dashes), and the NNLL QCD-like solution expanded 
in powers of 1/ ln(A|^//i^) given in (2^) (dots). The symbols and labels ai^e the same as in Fig. [I| For 
comparison, the LL (one-loop) solution (long dashes) is shown again. 



Mh are of key importance: If the Higgs self-interaction is non-perturbative at the scale Mh at 
which ( |1 ■ 1| ) is imposed, we do not expect perturbation theory to work in any other observable. 
Two-body Higgs decay rates are examples for such one-scale processes. They do not contain 
large logarithms, which need to be summed to all orders. Nevertheless as emphasized at the end of 



Sect, pi A| , RG methods can be used to judge the accuracy of perturbation theory: The dependence 
on the renormalization scale /i, at which the decay rate is evaluated, is of the neglected order in 
A. The stability of the perturbative result with respect to variations of fi must therefore increase 



order-by-order in A. This criterion to test perturbation theory has been first used in QCD in [ [15| ] 
and has become a standard method in QCD. Likewise the ratio of two results obtained in different 
renormalization schemes must approach unity with increasing order in A. This renormalization 
scheme dependence will be the second tool used in the analysis of the breakdown of perturbation 
theory. 



A. Higgs decay into gauge bosons 

At first we consider the decay rate of a Higgs into two gauge bosons. The 0(A") corrections 
to the decay rate H — * W^W^ are the same as the ones for the decay into two Z bosons. The 
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decay rates of the two channels only differ by an overall factor of proportionality. 
We write for the decay rate 



As 

A tree 



T{H ^ W^W-) = Ttree^^ (1 + AFscheme) • (4.1) 



Here Ttree is the Born approximation with 



tree 



^ 4M^ / 4M^ ^ 12M^\ A^_ 



The corrections stemming from the Higgs self- interaction have been calculated in the OMS scheme 
to two loops [|T6|]: 



1 + AFoMs = 1 + 2.8OA0MS + 62.15A^Ms + O (A^ms) + O ( ^ j . (4.3) 



In [|T6p it has been stressed that the two-loop OMS correction exceeds the one-loop OMS term for 
Mh > 930 GeV, although the one-loop correction is still small compared to the tree-level term. 
This indicates that either perturbation theory does not work for M^ > 930 GeV or the one-loop 
term is accidentally small in the OMS scheme. 

Using ( |33| ) we can express T in terms of the MS coupling: 

1 + AFhs = 1 - W + 5.88^ + Afg (l441n^ ^ + 297.0 In - 122.7^ 

+ 0(Ay + o(^). (4.4) 

Here /i is the renormalization scale at which the decay rate is evaluated. The renormalization point 
/io, at which the coupling is defined by in terms of Gp and Mh, is chosen as fio = Mh 
throughout this section. If we also take fi = Mh, we observe that the apparent convergence 
of AFj^ is worse than that of AFoms since both Aj^ > Aqms and the MS coefficients of the 
perturbation series are in magnitude larger than in the OMS scheme. The two-loop MS correction 
equals the one-loop MS term for Mh = 770 GeV. The corresponding OMS result yields Mh = 
930 GeV 

Let us now investigate whether one can refine these bounds on a perturbative Higgs mass by 
the examination of the renormalization scheme and scale dependence. 

We start our analysis by looking at the scheme dependence of the decay width. In the left 
plot of figure Fig. ^ we show the normalized scheme dependence (Fqms ~ ^ms) /^tree at one 
and two loops, using NLL and NNLL running couplings, respectively. The plot indicates that for 
Mh < 469 GeV the scheme dependence is reduced when going from one-loop to two-loop order. 
For larger Higgs masses, the scheme dependence increases, suggesting that perturbation theory 
is not working satisfactory or even fails. The scheme-dependence criterion indicates problems 
with perturbation theory in at least one of the two schemes considered for Mh > 469 GeV 
The criterion, however, is sensitive to the possibility of an accidental smallness of the one-loop 
correction. We remark here that the two-loop MS term in ( [4.4D is less than four times the square of 
the one-loop term, while in the OMS scheme, (|4.3D, this ratio almost equals eight. For Mh = 700 



16 



GeV the result for Toms is larger than the MS expression by 23%, and for Mh = 780 GeV the 
scheme dependence reaches unacceptable 52%. The simple criterion of comparing the magnitudes 
of the one-loop and two-loop contributions gives upper bounds on Mh which seem to be too large. 

Another criterion for the validity of perturbation theory is the order-by-order reduction of the 
scale dependence. The explicit yU-dependent logarithms in ( P| ) compensate the effect of the run- 
ning coupling to the order considered. To use this criterion in the OMS scheme as well, 
we need to introduce the RG logarithms into the OMS decay width (|4.3|). In the OMS scheme. 



the Callan-Symanzik equation describes the response of some Green's function to the scaling 
Pj — > fi/Mn ■ Pj of its external momenta. Its solution for the decay rate ( P3| ) reads 

1 + AToMS = 1 + AoMs(/i) (^-12 In + 2.80^ 

+ A^Ms(/^) (l44 In^ ^ + 88.8 In ^ + 62.15^ 

+ 0(aSms)+0(^). (4.5) 



Expanding the running coupling AomsIa*) in ( P~T| ) and ( P~5| ) in terms of Xoms{Mh), ( F3D yields 



( P3| ) up to the neglected order Aqms' of course. The OMS result expressed in terms of the running 
coupling ([4.5p can be used to examine the scale dependence of the decay width. 



We use the LL, NLL and NNLL expressions of the decay rates which consist of the Bom, 
one- and two-loop result in ( |4-.4D supplied with the LL, NLL or NNLL running coupling X^{fi), 
respectively. Then the renormalization scale fi is varied. The right plot in Fig. ^ shows that 
perturbation theory nicely works for the MS scheme if Mh = 400 GeV: The scale dependence 
diminishes order-by-order. A similar behavior is found for the OMS result using the OMS running 
coupling and the same value of Mh- Next we look for an upper bound for a perturbative Higgs 
mass. We investigate the scale dependence of the decay rate for values of Mh up to 800 GeV. 
Based on the size of the logarithmic terms in the higher orders of the perturbation series for the 



we have already advocated the range O.SMh < /i < 1.25Mh in Sect, pi A| . We may add a 
physical argument for this range as well: Suppose one decides to include the non-zero width of 
the Higgs into the analysis. Then the decay diagrams must be calculated with an off-shell Higgs 
boson with invariant mass ^/s, and the result is convoluted with a Breit-Wigner function. The 
decay rate would differ from ( |4-.4D by a function of s and Mh vanishing for s = Mh- The decay 
rate then involves two logarithms: ln(/i^/s) and ln(yU^/M^). The choice of fi could be with equal 
right /i = Mh, A* = or any scale in between. This suggests choosing the range for /i to be 
of the order of the total width. For the values of Higgs masses we are interested in, the width is 
between 0.2Mh and 0-3Mh [J^], so that the range O.SMh < fx < 1.25Mh is appropriate in both 
OMS and MS scheme. 

In Fig. ^ we have plotted the scale dependence of T(H W^W^)/ (Ftroe) vs. the physical 
Higgs mass in both OMS and MS scheme. The tree-level coupling is Atrcc = GfMh/V2- The 
scale dependence at a given order is represented by the smallest and the largest value of F /Ftree 
when ji is varied in the range Q.SMh < jJ^ < 1.25Mh- In the MS scheme the scale dependence 
correctly decreases when passing from LL to NLL to NNLL order if Mh < 742 GeV. For Mh = 
742 GeV the scale dependence in the LL and NLL order become equal and reach 36%. For 
Mh > 750 GeV and /i = 1.25Mh, the consistent solution of the NNLL running coupling is no 
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FIG. 4. The left plot shows the scheme dependence (Foms — rMs)/^ti'ee of the decay rate 
T{H — > W~^W~). For Mh > 469 GeV the two-loop scheme dependence (solid line) is larger than 
the one-loop scheme dependence (dashed). The right plot shows the dependence on the renormalization 
scale ^ of Fj^ for Mh = 400 GeV. The scale dependence diminishes sizably order-by-order indicating 
that perturbation theory works well for this value of Mh- 



longer defined in the MS scheme, indicating the ultimate breakdown. In the OMS scheme, the 
NLL and NNLL scale dependences become equal for Mh = 672 GeV, but are numerically small 
(8%). 

We conclude that perturbation theory for bosonic Higgs decays breaks down for Higgs masses 
of the order of 700 GeV. The scale-dependence criterion yields similar upper bounds on Mh in 
both schemes, although the absolute scale dependence is much smaller in the OMS scheme than 
in the MS scheme. Using running coupling solutions other than the consistent one (P3|), we obtain 
similar bounds. 

It should be noted that our scale-dependence criterion is not only sensitive to the coefficients 
of the different orders in AF, but also to the coefficients of the /5-functions which also enter the 
non-logarithmic terms of the uncalculated higher orders via diagrams connected with counterterms 
in (Oh. 



B. Higgs decay into two fermions 

Our next example is the fermionic decay width of the Higgs particle. At Born level it reads 

Here = (3) for lepton (quark) flavors. At tree level, this process is independent of the 
Higgs coupling A. Including radiative corrections, the fermionic decay rate receives corrections in 
powers of A. In the OMS scheme they are [P,[T7|] 
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FIG. 5. The scale dependence of T/Ttrcc for H W+W in the MS (left plot) and OMS (right 
plot) scheme. The extremal values of T/Ttrcc are plotted for O.SMh < ^ < 1.25Mh in the LL (dotted), 
NLL (dashed) and NNLL (solid) approximation. Except for the MS NNLL solution the extremal values 
correspond to the chosen bounds ^ = 0.8Mh and /i = 1.25Mh as indicated in the figures. The scale 
dependence diminishes order-by-order for Mh < 742 GeV in the MS scheme and for Mh < 672 GeV in 
the OMS scheme. The arrows indicate the equality of the LL and NLL scale uncertainty in the MS scheme 
(left plot), and the equality of NLL and NNLL scale uncertainty in the OMS scheme (right plot). 



1 + AFoMs {h ^ //) ^ 1 + 2.12A 



OMS 



32.66A0MS 



+ O (^OMs) • 



(4.7) 



Since the tree-level result of the fermionic Higgs decay, Eq. {\l.6\), is independent of the cou- 
pling A, we only need the one-loop relation between Aqms and A,^ to calculate the MS decay 
width up to two loops. Likewise, our scale variation criterion only involves LL and NLL running 
coupling in connection with the one-loop and two-loop results. As a result we can only compare 
the scale dependence of two instead of three orders. Since the LL one-loop result is identical 
in both schemes, the scheme dependence can only be compared at the NLL two-loop level. No 
bounds on Mh can be derived from the scheme-dependence criterion without knowing the three- 
loop corrections. This distinguishes the decay F (^H ff^ from the case F (i/ — > W^'*"!^") 
discussed in the preceeding section. 

Combining Eq. ( |3J| ) with the previous equation, we obtain the correction to the fermionic 



Higgs decay in MS quantities [U^JT^] : 

1 + AFms(^-//) ^1 + 2.12Ams- (51.03 + 25.411n(/iVMi)) A|s + 0(A|s; 
We also give the scale dependence of the OMS result: 



1 + AF 



OMS 



1 + 2.12AoMs(yu) - 32.66 + 25.41 In 



M'h, 



^OMS 



(/i) + o (a 



3 

OMS 



(4.8) 



(4.9) 



In Fig. ^ we show the scale dependence of the decay width expressed in terms of the running 
coupling in both MS and OMS scheme at one and two loops as a function of Mh- At each order 
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FIG. 6. The resummed /^-dependence of the normahzed fermionic Higgs decay width. Only the 
//-dependence entering through the Higgs couphng is analyzed, i.e. no running Yukawa or QCD couplings 
are considered. Using the MS renormalization scheme (left plot), the /^-dependence of the 2-loop result is 
as big as the //-dependence of the 1-loop result for Mh = 513 GeV as indicated by arrows. In the OMS, 
this equality occurs for Mu = 776 GeV (right plot). The scale dependence is nevertheless small at these 
critical Higgs mass values. 



the three curves refer to // = O.SMji^, X.^Mu, and 1.25Mh- We find that the fi dependence of 
the two-loop result is larger than the /i dependence of the one-loop result if is larger than 513 
GeV in the MS scheme. In the OMS scheme the corresponding bound is Mh = 776 GeV. These 
results, however, are not as valuable as those found in the previous section, since they are only 
founded on the comparsion of two orders rather than three. The low value of Mh = 513 GeV in 
the MS scheme may be accidental. The scale dependence is very weak, less than a few percent for 
Higgs masses up to 750 GeV This is due to the fact that the tree-level result does not depend on A. 
The scheme dependence at NLL (two-loop) is marginal for the same reason. We presently have no 
information on the NNLL behaviour of the fermionic decay. Looking at Fig. ^ we conclude that 
the upper bound on a perturbative Higgs mass is in agreement with our findings in the case of the 
bosonic Higgs decay. 



V. SCHEME AND SCALE DEPENDENCE OF CROSS SECTIONS 



The previous section considered one-scale processes. The only logarithms appearing in the 
RG improved results are ln(yU^/M^). Testing the perturbativity of Higgs sector is much more 
stringent when considering high-energy two-scale processes, because they involve the running 
coupling at a high scale. Typical two-scale processes which depend on the coupling A are 2 2 
scattering processes involving longitudinally polarized gauge bosons and the Higgs boson. In 
the limit s, Mfj ^ the longitudinally polarized gauge bosons can be replaced by massless 
Goldstone bosons. In this limit all couplings except for the trilinear and quartic Higgs coupling 
are subleading and can be neglected. Also assuming s 3> Mfj, the high-energy amplitudes and 
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cross sections have been calculated to two loops by [|8|J20|]. The high-energy limit is approached if 
y/s > 2Mh [pot]. Its error is less than a few percent if y/s > 5Mh- For most of our analysis we 
will consider y/s = 2 TeV. 

A typical example is the process W^W£ ZlZl- The high-energy cross section only 
depends on i/i, Mh, and A. The OMS cross section is p^ ] 
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where the OMS renormalization fixes /xq = Mh such that the OMS running coupling has the 
boundary condition Aoms(/^ = /^o) = Mfj/{2v'^). The scale /i must be chosen of the order to 
resum the large logarithms to all orders. For a complete resummation of all logarithms up to three 
loops, one uses the NNLL running coupling and takes ^ = y/s. The yU-independent logarithm is 
related to the field anomalous dimension and needs no resummation because of the smallness of 
its coefficient. At one loop, the anomalous dimensions of the Higgs and gauge bosons are zero, 
and at two loops they are numerically unimportant [jSj]. 
The MS resuh is 
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We choose the renormalization point to be fiQ = such that the the MS running coupling is 
fixed at /i = /io by (PT^f). 

The RG structure of the cross section is similar to the one of the decay width T(H ^ W^W~) 
since the tree-level result also depends on A. There are two important differences: the tree- 
level cross section is proportional to A squared, and the running coupling resums terms of or- 
der ln(s/M^) which can lead to a significant increase of the running coupling compared to the 
tree-level coupling. 



A. Perturbativity at collider energies 

Choosing ^ = yfs = 2000 GeV, we show the scheme dependence of the cross section 
^L^L ~^ ZlZl in Fig. |^. For Higgs masses larger than 400 GeV, the scheme dependence 
is larger than 40% and actually increases when going from NLL to NNLL cross section. Choosing 
fi = y/s = 1000 GeV, a value more realistic for lyiy-scattering at future colliders, we find the 
critical mass value to be M^ = 436 GeV. For such values of M^ and y/s the high-energy approx- 
imation yields about 70% of the exact A-dependence of the cross section [|20|]. Increasing y/s, the 
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breakdown of perturbation theory as seen in the scheme-dependence criterion happens for smaller 
and smaller values of Mh- We will come back to this later. 

In addition to the scheme dependence, we also evaluate the scale-dependence criterion. Since 
the renormalization group is used to resum ln(s) terms, the scale /i is varied around ^/s rather 
than Mh- The arguments of Sect. 0, which are based on the size of the beta function coefficients, 
suggest the range < fi < 1.25y^. Choosing = 2000 GeV, the result of this variation 

is shown in Fig. ^ for both MS and OMS scheme. In the MS scheme, we observe a nice order- 
by-order reduction of the scale dependence if Mh is less than 357 GeV. For larger values of Mh, 
the NLL scale dependence exceeds the LL one, and for Mh > 368 GeV we also find the NNLL 
scale dependence to be larger than the NLL one. For the OMS result we find only one crossing 
point. There the NNLL scale dependence dominates over the NLL one if Mh > 410 GeV. Taking 
also the scheme dependence into account we conclude that perturbative results for longitudinal 
gauge boson scattering at = 2000 GeV cease to be meaningful if Mh is of the order 400 GeV 
or above. For a value of ^/s = 1000 GeV we find limits of 0(450 GeV), but here the additional 
low-energy contributions — though not dominant — may change the limit. 

The upper bounds on a perturbative Higgs mass from scattering processes are more stringent 
than the results found in Higgs decays. The reasons are two-fold. First, the one-loop and two-loop 
coefficients of the perturbative corrections to the cross section are an order of magnitude larger than 
the corresponding coefficients of the decay processes. Second, the running coupling evaluated at 
a scale y/s > Mh is numerically larger than the coupling involved in the decay widths. 

It is also interesting to note that the bounds derived here are similar to the results found in [ ]21| ] 
which are based on perturbative violation of unitarity at two loops, a rather different criterion for 
the breakdown of perturbation theory. Compared to the nonperturbative lattice results of Mh < 
710 ± 60 GeV []22|], our perturbative criteria require smaller Higgs masses. 



B. Perturbativity at large embedding scales 

In Fig. ^ we show the upper bound for a perturbative Higgs mass as a function of some embed- 
ding scale A. Since the Standard Model Higgs sector is — depending on the value of Mh — not 
defined above a certain energy scale, it is a common procedure to introduce an embedding scale 
above which the physical interactions are to be described by a more complete theory. Requiring 
that the SM Higgs sector is still perturbative at such an embedding scale, it is possible to give up- 
per bounds on the Higgs mass as a function of the embedding scale. Here we require the process 
Wi^Wl — > ZlZl to be perturbative for energies ^/s < Aembed and apply our criteria for scheme 
and scale dependence to calculate the upper bound on Mh- The result is shown in Fig. ^. At 2 TeV, 
the different bounds on Mh correspond to the values derived from Figs. [7] and For increasing 
embedding scale, the upper bound on Mh decreases, and the three different criteria give converg- 
ing bounds. At Acmbed = 10^^ GeV, the upper bound is 150 ± 3 GeV. For such low Higgs masses, 
however, one expects the top-quark Yukawa coupling to have an influence on the RG evolution of 
the Higgs coupling. As a matter of fact, the SM beta function of the Higgs coupling can become 
negative if the Higgs particle is too light, invalidating our analysis. Taking rrit = 175 GeV, this 



is expected to happen for Mh < 135 GeV [p30, which is lower than the values considered by us. 



Using the results of Lindner [E^, we expect the change of our Higgs mass bounds due to a top- 



quark mass of 175 GeV to be less than 10% for the whole range of energy-scales considered. For 
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FIG. 7. The scheme dependence of the cross section W^Wi^ —>■ ZlZl at ^/s = 2 TeV. For 
Mh > 400 GeV the two-loop scheme dependence (solid line) is larger than the one-loop scheme de- 
pendence (dashed) and exceeds 40%. 



comparison, we also display the values of Mh which would lead to a one-loop Landau singularity 
at Acmbcd- Calculating the Landau pole, we again neglected Yukawa and gauge couplings. 

It is worth noting that our criteria for an upper Higgs mass lead to values of the running 
coupling at the embedding scale which are not very large. For values of ^/s in the TeV range, 
the maximal allowed value of the running coupling is less than 2.2, and at GUT energies we find 
the maximal value to be less than 1.6. These small values are surprising: The beta function and 
the solutions for the running coupling show excellent convergence for such small values. As a 
matter of fact, the usual criteria for the breakdown of perturbation theory assume that the running 
coupling becomes large. Yet this method is unsatisfactory: The consideration of the one-loop 
Landau pole as in cannot be extended to higher orders. Instead the authors of ^ use the 
criterion of (in our notation) A(Aembed) < tt^. Here the choice of the numerical bound involves 
some arbitrariness. The plots in Fig. [l| and Fig. ^ show that vr^ is clearly chosen too large, because 
the different solutions differ substantially for A = vr^. The second arbitrariness of this methods 
resides in the fact that the breakdown of perturbation theory is judged from a scheme dependent 
quantity: A(Aembed) depends on the renormalization scheme through the radiative corrections to 
the matching condition ( |1 . 1| ) and, most important, through the coefficients n > 2, of the P- 
function. In contrast physical observables are scheme independent up to the calculated order. They, 
however, seem to become non-perturbative for much smaller values of A. For comparison, the 
unitarity arguments used in [ PH ] yield an upper bound on the running coupling of 2.3, independent 
of any choice of ^/s. 
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FIG. 8. The resummed //-dependence of the normahzed W^W^ ZlZl cross section at -s/s = 2 
TeV. Only the /x-dependence entering through the Higgs coupling is analyzed, i.e. no running Yukawa or 
QCD couphngs are considered. Using the MS renormalization scheme, only a Higgs mass of less than 357 
GeV guarantees an order-by-order reduction of the scale dependence. In the OMS scheme, the NLL and 
NNLL are equal if Mh = 410 GeV. The scale dependence is almost 50% at this value of Mh. 



VI. CONCLUSIONS 

Our analysis of various physical processes of the Higgs sector indicates that the breakdown 
of perturbation theory cannot be judged purely on grounds of a large Higgs (running) coupling. 
The breakdown of perturbation theory in both Higgs decays (Aqms = CpM^/y/^ < 4) and 
scattering processes (X{^/s) < 2.2) occurs for relatively small values of the Higgs coupling. The 
usual criterion — the breakdown of the perturbative behaviour of the beta function and the running 
coupling — yields upper bounds for a perturbative Higgs mass which are too large. 

Applying our criteria of scheme-dependence and scale-dependence to Higgs decays, we find a 
satisfactory perturbative behaviour of the decay widths if Mh < 0(700 GeV). In the case of 2 — > 2 
scattering processes, one needs to specify the energy scale at which the process is to take place. 
Choosing ^/s to be a couple of TeV, the Higgs mass has to be less than 0(400 GeV) to guarantee a 
good perturbative behaviour of the cross section for W^Wj^ ZlZ^ scattering. Requiring such 
processes to be perturbative at y/s — 10^^ GeV, the Higgs mass has to be less than 0(150 GeV). 
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FIG. 9. Upper bound on Mh derived using our scheme-dependence and scale-dependence criteria. 
We require the SM process W^W£ — > Z^Zl to remain perturbative for cms-energies up to ^/s = Agmbed- 
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